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Abstract 

The production of top quarks in single mode at future ep colliders is studied, 
the attention being mainly focused to the case of the proposed LEP©LHC 
collider. We are motivated to reanalyse such a process following the dis- 
covery of the top quark at Fermilab. Thanks to the measurement of its 
mass one is now able to establish more accurately the relevance of single 
top production for itself and for many other processes to which it may act 
as a background. In addition, the recent improvement of our knowledge of 
the quark and gluon dynamics inside the proton now allows one to pin down 
the dependence of single top production on the partonic structure functions. 
Both the leptonic and hadronic decay channels of the top quark are studied 
and compared to the yield of the corresponding irreducible background in 
presence of 6-tagging. 
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1. Introduction 



Now that the Fermilab experiments have clearly assessed the existence of the heaviest 
quark of the Standard Model (SM) and given a rather accurate measurement of its 
mass 0, many of the theoretical calculations carried out in the previous years need to 
be updated to the current value of this fundamental parameter. In this paper, we turn 
our attention to the case of top production in single mode at future electron(positron)- 
proton colliders. As a further motivation for our revision we put forward the fact that 
a huge amount of data improving our knowledge of the parton distribution functions 
(PDFs) has been produced in the years following the early studies of top phenomenology 
at ep machines (see, e.g., Ref. 0]), along with more detailed treatments of the dynamics 
of heavy quarks inside the proton. Therefore, the error associated with the partonic 
behaviour in the initial state should be at present significantly smaller than in the 
past. Finally, the reduction of the theoretical uncertainty in single top production also 
implies that other effects, such as those due to the irreducible backgrounds, need now 
to be incorporated in more detailed phenomenological analyses. We calculate these 
effects here for the first time, in both the hadronic and leptonic top decay channels. 

In order to illustrate the particular relevance of single top processes in 
electron(positron)-proton annihilations we remind the readers of the motivations for 
higher energy ep experiments. 

Firstly, such colliders will be an obvious and unrivaled testing ground for QCD at 
very low Bjorken a; [Q, in exploring the structures of both the proton and the photon |^ 
at the TeV scale, taking over the presently running HERA machine 0. In connection 
with this point, we will show that the single top quark process discussed here can be 
be useful in understanding the phenomenology of the PDF of the bottom quark. 

Secondly, and particularly in the case of the proposed LEP©LHC collider P|, they 
will be able to search for the Higgs boson of the Standard Model [0 (or the lightest 
neutral Higgs boson of the Minimal Supersymmetric Standard Model, MSSM) in the 
intermediate mass range 90 GeV ~ ~ 130 GeV P, P[, in the case it may not be 



accessible at the LHC nor the existing colliders (see Ref. |T0| for discussions on this 
point). As was discussed in the jet background due to the top quark will be large 
if only a single 6-tagging is implemented in identifying the Higgs boson decay in the 
most favoured channel ep —>■ VeW^W~ X —>■ z/e0X — >■ UgbbX 

Thirdly, the role of an ep machine will be complementary to those of e~^e~ (i.e., 
NLC) and pp (i.e., LHC) colliders in the search for New Physics such as leptoquarks, 
excited leptons, low mass sleptons, doubly charged Higgs bosons, and new vector bosons 
(see @] and the references therein). Many of these processes have neutral current-type 
interactions of the form eg — > eg, and so the single top quark process, when the top 
quark decays leptonically to a bottom quark, a positron and an electron neutrino, is 
a potentially dangerous background which should be included in experimental simula- 
tions. As for the possibility of exotic top quark decays, the study of the Supersymmetric 
two-body decay modes is a straightforward extension of this project and will be 
carried out elsewhere [|12]. 

Although the physics potential of a higher energy ep machine is suppressed compared 
to a pp one by the reduced centre-of-mass (CM) energy and luminosity, we stress its 
allure in the suppression of the initial state QCD noise, which allows for a cleaner 



environment to study the physics of the TeV scale, possibly before an NLC will be in 
operation We also mention that the physics of ep colliders, in conjunction with 



the discussed possibility of their running in the jp mode |T|], has been recently under 
renewed and active discussion ITI 



The production of top quarks at future ep colliders has been studied in the 
context of top quark searches at LEP1©LHC during the 1990 Aachen Workshop. 
A detailed study was presented in the corresponding Proceedings fl^. There, the two 
following channels were investigated 1171 : 



CC: W^ + g^ tb, (1) 

NC: j,Z + g^tt, (2) 

via charged (CC) and neutral current (NC) scatterings of an off-shell gauge boson 
against a gluon, the former being produced via bremsstrahlung off the incoming elec- 
tron(positron) and the latter being extracted from the proton. In general, the CC 
channel dominates over the NC one, due to the larger phase space available. For 
rrit = 175 GeV p|, the suppression is more than one order of magnitude at the TeV 
scale (Indeed, this is the reason why we will concentrate on W^g fusion only.) 



In Ref. |jT6|, also a detailed signal-to-background analysis was carried out, in both the 
(semi)leptonic and hadronic top decay channels. 

The Feynman diagrams describing reaction (0) induced, e.g., by positron beams, 
can be found in Figure la, where the top is considered on-shell. As the bottom mass is 
small compared to rrit, the dominant contribution to the total cross section comes from 
diagram 1 of Figure la, when the final b quark is coUinear with the incoming gluon. 
The collinear divergences are however regulated by the finite value of the bottom mass 
and manifest themselves by means of contributions of the form L = as{n'^) log(/x^/m^), 
with /i^ ~ s, being s the CM energy at the 'partonic' level e^b and the strong 
coupling constant. 

Such logarithms are rather large, thus terms of the form a"(/i^) log" (/i^/m^)/n! have 



to be resummed to all orders in perturbation theory |]18| in order to compute the cross 
section rehably. This can be done by introducing a b parton distribution, fb{x,iJ?'), in 
terms of the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) splitting function 

P,,{z) = \[z^ + {l-z^)] (3) 

of a gluon into bb pairs {z being the fractional energy carried away by the (anti)quark). 
In fact, bottoms are not valence quarks, rather they materialise once the energy scale 
H of the evolution reaches their production 'threshold' at a given value Hh ~ "^b. The 
function Pbg{z) is indeed the 'coefficient function' of the logarithmically enhanced term. 
The b structure function then evolves with /i according to the DGLAP equations, from 
an initial condition of the sort, e.g., /f,(x,/i^) = if /x^ < fif. 

It follows then that single top production and decay via process (|I]) can be conve- 
niently studied by computing the transition amplitude squared for the reaction (e.g., 
assuming incoming positron beams) 

e+b ^ Uet ^ i^ebW^ Uebff, (4) 



where / represents a lepton/neutrino or a u,d,s and c quark (produced in the top 
decay), appropriately convoluted with a b distribution function evaluated at the adopted 
scale /i^. We exploit here this approach. 

In our opinion, such a procedure (in which the parton is b) is more appropriate 
than the one exploited in Ref. |T^ (in which the parton is g), especially at high en- 



ergies. In fact, we have explicitly verified that for the values of considered here 
the dominant contribution to process (|I]) comes from configurations in which the b is 
collinear with the incoming gluon. Since it is exactly such emission that is summed to 
all orders in perturbation theory in leading and next-to-leading logarithmic accuracy 
inside the b structure function when ^^(/i^) log(/i^/m^) ~ 0(1) (and this is clearly the 
case in our context, e.g., when /z^ = s ~ m'f), our approach will give a more accurate 
answer. However, for comparison, we will also show in the present paper several rates 
as produced by the process induced hj g bb splitting^ 

In this paper we study the single top quark production via e^b fusion at various 
energies, together with all tree-level irreducible background processes as shown in Fig- 
ures lb, Ic and Id. Figure lb corresponds to the case of the leptonic decays of the 
boson in the signal process, 

e+6 ^ Uebti^i, (5) 
where i = e, fi, r, whereas Figure Ic refers to hadronic decays of the boson, 

e+b ^ uebi+jf, (6) 

where jj' represents a pair of light quark jets ud or cs. To these must be added the 
case of the gluon mediated background of Figure Id, 

e+b ^ uM'-jf, (7) 

where jj' again represents a pair of light quark jets. 

In addition to these three, if charge measurements of the bottom quarks prove 
impractical or impossible, we will have background from processes of the forms: 

e+b ^ Uebi+i^i, (8) 

e+b ^ uM^jf. (9) 

A single 6-tagging capability is assumed throughout, all results being linearly pro- 
portional to its efficiency. The case i = t assumes that jets coming from the tau and 
the quarks will easily be distinguishable. We perform all calculations for the case of 



^Note that the complete next-to-leading (NLO) corrections to W+g/b fusion involving the collinear 
logarithms as well as the large angle emission and the loop diagrams have been recently presented in 
the MS renormalisation scheme jl^ . Such results contradict earlier ones based on the DIS factorisation 
scheme For the case of ep collisions at HERA they amount to approximately 2% of the result 

obtained by means of the b structure function approach, and they are rather insensitive to collider 
energies in the TeV range. Therefore we expect them to be well under control also at the proposed 
LEPqLHC, so for the time being we do not include them in our calculation. Another reason for 
doing so is that we will also be concerned with the interplay between the single top signal and the 
non-resonant irreducible background, which is here computed at lowest order. 



e~^p colliders, although the e~p case is precisely identical since the calculations involve 
no valence quarks and are therefore invariant under the exchange ^ e~ . 

The plan of this paper is as follows. In Section 2, we describe the methods we 
adopted in the calculations of the signal and background processes. In Section 3, we 
present and discuss our results. Section 4 is a brief summary. 

2. Calculation 

The tree- level Feynman diagrams that one needs in order to compute processes (|^), (j^) 
and are given in Figures lb, c and d, respectively. For reaction (|^) we show the 
diagrams for the case i = e, which is the most complicated. When i = fi or t only 
ten out of the twenty-one diagrams in Figure lb contribute. For processes the 
number of diagrams is independent of the flavour j'Q. 

The single top quark signals @ are produced by diagrams 11 in Figure lb and 4 in 
Figure Ic, for leptonic and hadronic decays, respectively. The remaining diagrams 
in Figures lb and Ic represent the 'irreducible' background to single-top production and 
decay. Reaction (|^) does not contribute to the signal at all, but only to the background. 

Graphs in Figures lb, c and d refer to the case of e^b fusion, i.e., to the scattering 
of a positron and a bottom quark, the latter being extracted from the incoming proton 
beam. As mentioned earlier on, we have treated the bottom quark as a constituent 
of the proton with the appropriate momentum fraction distribution /^(x, /i^), as given 
by our partonic structure functions. It can be noted that the bottom antiquark is also 
present inside the proton with an equal probability. When calculating rates for single- 
top production at e~^p colliders, diagrams initiated by bottom antiquarks must also be 
considered. However, as long as the deep inelastic scattering of the proton takes place 
against a positron, such graphs do not produce a resonant top quark. The topologies of 
these bottom antiquark initiated graphs are easily deducible from those in Figures lb 
and c. From the point of view of top quark studies, these act as additional backgrounds. 
Their production rates will be different from the case of e^b fusion if the CM energy 
at partonic level (i.e., spans the top quark production threshold. In contrast, the 
cross sections due to e^b initiated diagrams and proceeding via QCD interactions are 
identical to the yields of reaction (0) and the actual graphs are the same as those in 
Figure Id, apart from the trivial operation of reversing the bottom quark line^ 

The possibility of the '26' charged current processes e'^q Uebbq (where q is d, u, 
s or c) being mistagged as a single b event and acting as background to the hadronic 
channel can not be neglected, even when the 6-tagging efficiency is high. If the latter is 
denoted by e^, then the probability of misidentification is given by 26^(1 — e^), assuming 
no correlations between the two 6-taggings. Thus the suppression of the 2b background 

■^We will refer to process (||) as the 'leptonic' channel, and to processes (^) and (|^) as the 'hadronic' 
channels. For the latter cases, we will further distinguish between 'electroweak' (EW) and 'strong' 
(QCD) production, respectively. 

^Note that for the case of e~p scattering things work in a complementary way, the resonant top 
antiquarks being produced by incoming b partons. Indeed, as bottom (anti)quarks are produced inside 
the nucleon via a g bb splitting (that is, they are sea partons), no differences occur in the deep- 
inelastic dynamics of the above processes if antiproton beams are considered. Although we study 
positron-proton colliders here, our discussions are transposable to all the other cases. 



with respect to the single b events is 2(1 — e?,). As our investigations concern mainly the 
top quark signal process @, the complete analysis of single bottom quark processes 
being outside the scope of our present study, we content ourselves with an estimate 
of the degree to which this additional 'irreducible' background could affect the top 
quark and boson mass reconstruction (the 2b background does not contribute to 
the leptonic case). Our explicit calculation, using the methods explained below, shows 
that after cuts in the reconstructed top quark and boson masses are introduced, 
the cross section of the 2b process is of the same order as that of reaction (|^) (differing 
only by 10% at the LEP©LHC energies), thus being quite small in the end (see Section 
3.). 

To calculate the squared amplitudes for processes (^-(^ we have used the FORTRAN 
packages MadGraph [pl|] and HELAS p2|. The codes produced have been carefully 



checked for gauge and BRS invar iance at the amplitude squared level. The multi- 
dimensional integrations over the phase spaces have been performed numerically using 
the Monte Carlo routine VEGAS p4|, after folding the partonic differential cross sec- 



tions with the appropriate quark densities. The programs that we have produced have 
been run on a DEC 3000 Model 300 alpha-station, on which the evaluation of, e.g., 
10^ events took some 14 minutes of charged CPU time to produce a cross section at 
the level of percent accuracy in the case of process (|^) for the sum of the two contri- 
butions i = e and fj, (the latter being equivalent to the case £ = r): that is, for the 
channel involving the largest number of diagrams and the most complicated resonance 
structure. 

All the codes implemented are available from the authors upon request. To allow 
for a prompt evaluation of single-top rates at any energy and for any choice of selection 
cuts, we have also calculated the amplitude squared of process (^) analytically, including 
top width effects. In the leptonic case, and assuming all lepton and neutrino masses to 
be zero, it reads as follows^: 



\Me+b^p,b'£+ue\ = 2(47raem/ sin Ow) \Pw*\ \Pw\ \Pt\ Pb ■ Pp^-Pe ■ PePt + '^Pe ■ PtPe ■ Pt) 
with 

Pw* = 1/{pIv* - M^±), Pw = l/{Pw - + tMw±Tw±), Pt = l/{Pt - m^t + ^^tTt) 

PW*=Pe-Pu,, PW=P£+Pue, Pt = Pb + PW ■ 



In the hadronic case, again assuming zero light quark masses, the above formula needs 
to be multiplied by the colour factor 3, and and z/^ replaced by d{s) and m(c) 
respectively. 

As the default set of PDFs we used the NLO set MRS (A) |2^. However, as one of the 
motivations of this study is to investigate the dependence of process (^) on the evolution 
of the structure functions of bottom quarks inside the proton, we have produced our 

^The analytic expression for process (|l|), also involving the decay currents, can be found in Ref. |p5|| . 
We have checked our ME for the gluon induced process against that given in Ref. |^ in the appropriate 
configuration (i.e., for pp collisions) and found perfect agreement. 



results for other 23 recent NLO PDFs which give excellent fits to a wide range of 
deep inelastic scattering data and to others from different hard scattering processes 
(see the original references for details). These are the packages MRS(A', G, J, J', Rl, 
R2, R3, R4), MRS(105, 110, 115, 120, 125, 130), MRRS(1,2,3) and CTEQ(2M, 2MS, 
2MF, 2ML, 3M, 4M) g^, m E|, |1], 0, g. Note that in each case the 



appropriate value of Aqcd = was used. In particular, for the MRS (A) set we 

adopted A^ = 230 MeV. Unless otherwise stated, the QCD strong coupling constant 
as entering explicitly in the production cross section of process and implicitly in the 
PDFs was in general evaluated at two-loop order at the scale fi = a/I. The same choice 
has been made for the scale of the structure functions. The spread of the results as 
obtained from the different packages with respect to our default MRS (A) value (rather 
than the errors of the numerical integrations) can be taken as a possible estimate of 
the uncertainty of our predictions throughout the all papeiQ. 

The bottom quark sea distributions are not measured by experiment, but are ob- 
tained from the gluon distributions splitting into bb pairs by using the DGLAP evolution 



equations [Q. Therefore, the b structure functions are different from the light quark 
distributions, which do need to be measured as they involve non-perturbative QCD, 
for which a consistent theoretical framework does not exist. In contrast, the PDFs 
of b quarks evolve at energies of the order of the fermion mass or larger, so that 
their dynamics can be calculated by using the well assessed instruments of perturbative 
QCD. That is, given the PDFs of the gluon and of the light quarks, those of the b are 
precisely determined, as they do not contain any free parameters (apart from rub, of 
course) . 

We think that by the time that ep colliders at the TeV scale will begin to be 
operative, the uncertainty on the gluon distributions at medium and small x may be 
expected to be significantly smaller than at present, principally due to forthcoming 
improved measurements of the small x deep inelastic structure functions at HERA, 
and of large pt jet and prompt photon production at the pp (Di-)Tevatron at Fermilab 
and the pp LHC at CERN (the latter being scheduled to start running around 2005)|]. 
Therefore, detailed studies of single top events produced in electron-proton collisions 
will allow one to constrain the error related to the dynamics of the g ^ bb splitting in 
the DGLAP evolution. In fact, we expect the experimental information on b structure 
functions as collected at the end of the HERA, Fermilab and LHC epoch to be rather 
poor, if not inexistent. On the one hand, at the CM energy typical of the ep accelerator 
now running at DESY (-y/igp = 314 GeV) the content of b quark inside the scattered 
hadron is very much suppressed per se On the other hand, at both the Tevatron 
{\fSpp = 1.8 — 2 TeV) and LHC (^/ip^ = 10 — 14 TeV) the study of b induced processes 



inevitably proceeds through either the production of top quarks in single mode |^], 
whose signatures suffer from a huge background due to tt production via qq and gg 
fusion, or via pure QCD interactions, biased by a large amount of light quark and gluon 
jet noise. These two problems can in principle be solved by future ep colliders. Firstly, 
they will be operating at the TeV scale thus allowing for a very much enhanced content 



^We have verified that differences in the results similar to those obtained in case of process (^) also 
occur for the complete tree- level reactions (^-(0)- 

^In fact, the typical x values probed via process e.g., at the LEP®LHC, are of the order rn^/s 
or more, that is above 10~^, where the gluon density is already well known at present. 



of initial b quarks, which can be probed in the 'kinematically' more defined context 
of a DIS process of an electron (positron) against a proton. Secondly, as discussed 
previously, the single top mode via e^p collisions has a much larger cross section than 
top-antitop production induced by •jg and/or Zg fusion ||T6[| . 

As for the dependence of process (|) that one expects on the different b structure 
functions, it is worth reminding the reader some peculiar features of the sets considered 
here. For starting, whereas in the earlier MRS sets (excluding MRRS(1,2,3)) the bottom 
density is set to zero below threshold (i.e., fh{x,fi'^) = for /i^ < fil) and for /x^ > fif 
the bottom distribution is evolved assuming a massless quark, nib = 0, the most recent 
ones (i.e., MRRS(1,2,3) [^) implement a formulation which allows heavy quarks mass 
effects to be explicitly incorporated in both the coefficient and splitting functions in the 
parton evolution equations and the two regions /i^ ~ /i^ and /i^ ^ /i^ can be treated 
consistently. Indeed, an alternative approach is to treat the bottom as a massless parton 
above /i^ = fif that is, the mass effects are neglected in the splitting functions. 



although they are included in the coefficient function at NLO. The difference between 
the two techniques clearly resides in the fact that the neglected ml effects give NLO 
contributions during the DGLAP evolution^. (The approach of Ref. |^ has however 



been criticised in Ref. |^ , where it was made the point that the presence or absence of 
heavy quark masses in the DGLAP evolution kernels has no effects on the measurable 
cross section.) Furthermore, also the choice of the threshold Hb can vary, being in some 
instances set at /i^ = ml (see, e.g., Ref. |2^) and in other cases at ^l = 4ml (^^e, e.g., 
Ref. m). 

For consistency with the parameter values adopted in the set MRRS(1,2,3) we have 
used as default for the mass of the bottom quark mb = 4.3 GeV [Q. (Note that the 
charm mass in the three above packages has been set equal to 1.35, 1.50 and 1.2 
GeV, respectively: this is however a 'dummy' value in the production process (^).) For 
the top mass we have taken (unless otherwise stated) mt = 175 GeV whereas for 
the width Ft we have used the tree-level expression [^. Leptons and u,d,s (and c 
as well) quarks were considered as massless in processes (^-(0). For simplicity, we 
set the Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix element of the top-bottom 
coupling equal to one, the Standard Model prediction at the 90% confidence level 
^ being 0.9989 < \Vtb\ < 0.9993. For the gauge boson masses and widths we used 
Mz = 91.19 GeV, Tz = 2.50 GeV, Mw± = 80.23 GeV and Tw± = 2.08 GeV. The 
electromagnetic coupling constant and the weak mixing angle are aem = 1/128 and 
sin^ 9w = 0.2320, respectively. 

The Higgs boson of the Standard Model enters directly in the diagrams of Figures 
lb (graphs 3 and 16) and c (graph 7), when the bottom quark mass is retained in the 
fermion-fermion-scalar vertex. As default value for the scalar mass we used Mh = 150 
GeV, according to the best fit as obtained from the analysis of the LEP and SLC high 



precision EW data: i.e., Mh = 149t^f GeV Q. However, since the constraints on the 



Higgs mass are rather weak (a lower bound of 66 GeV from direct searches and a 
95% confidence level upper limit of 550 GeV from the data mentioned |^) we studied 
the Mh dependence of the EW contributions in processes (§|)-(§) and (§|)-(y), and found 
it negligible (note that the Higgs boson is always produced via non-resonant channels 

^Note that a third approach exists in Htcrature which does not treat bottoms as partons. For 
example, b quarks are not present in the GRV sets of PDFs [Eol . 



in those reactions). This is also true for the 2b process, after the implementation of the 
selection cuts (see below). 

Finally, as total CM energy \/Sep colliding positron-proton beams we have 

adopted values in the range between 300 GeV (i.e., around the HERA value) and 2 
TeV. However, we focused our attention mainly to the case of a possible LEP2©LHC 
accelerator, using a 100 GeV positron beam from LEP2 and a 7 TeV proton one from 
the LHC, yielding the value y/s ~ 1.7 TeV in the CM frame of the colliding particles. 



3. Results 

As emphasised in Section 1, we generate the single top quark in the final state by means 
of the matrix element for e~^b fusion (^) convoluted with b structure functions rather 
than producing the initial b quark via an exact g ^ bb splitting folded with a gluon 
density. However, to investigate the differences between the two procedures, we show in 
Figure 2 the total cross section of the signal process (H) plotted against the CM energy 
of the ep system along with the yield of reaction (|I|) (the latter including top decays 
and finite width effects on the same footing as the former )0. 

Care must be taken when comparing processes (|l]) and (^ with respect to each 
other. In fact, one should recall that the corresponding rates are strongly dependent 
on the (factorisation) scale fi. In general, the W^g fusion cross section decreases sharply 
as the scale increases, whereas that of e~^b events goes up mildly as fi gets larger (see 
Ref. for a dedicated study in the case of pp collisions at the Tevatron). Although 



at LO there is no privileged choice for /i, Ref. |T9[ has shown that the most appropriate 
scale at the exact NLO (when both processes (0) and (^ need to be calculated) in the b 
distribution function is /i^ ~ + m^, where = —q^ {q being the four-momentum of 
the incoming virtual boson). Therefore, we have adopted this value in producing 
Figure 2 (also as argument of the strong coupling constant), whereas in all other cases 
we will maintain the LO 'running' choice /z = y/I. This has been done for two reasons. 
First, we have verified that for /i ~ Vlmin ~ fnt the rates of process (^) are rather 
stable, showing variations below 6-7%. Second, this choice of the scale allows one 
to consistently incorporate the non-resonant diagrams along with the top ones when 
calculating the cross sections of the complete processes (|)-(0). 

From Figure 2, it is clear that, apart from the different normalisation, the threshold 
behaviour in processes (|lD and (^ is substantially similar as a function of the total 
CM energy. Though, the ratio between the two series of curves is approximately 4.7- 
5.3 at 300 GeV and it decreases with increasing energy, stabilising at 1.3 TeV or so 
around 1.7-1.8. We trace back the behaviour at large energies as due to the fact that 
the term ^^(/x^) log(/x^/m^) becomes constant because of large logarithms cancelling 
each other (^^(/i^) is in fact proportional to 1 / log{fi'^ / A'^^^^)) . In contrast, at smaller 
energies (well below the TeV scale) this is no longer the case and, in addition, graph 2 
of Figure la becomes strongly suppressed, thus explaining the increase of the observed 

^Note that in order to obtain a gauge invariant cross section for process in presence of a finite 
value of Tt we need to consider a set of three diagrams. That is, the two with resonant top production 
(i.e., those in Figure la with the additional decay t bW"^ — > bff) and a third one in which the 
ff current is attached to the off-shell fermion propagator in one of the graphs of Figure la. 



ratio. The value of the latter between the two cross sections when ^/iep -1-3 TeV can 
be understood in terms of the large logarithms entering in the resummation of the b 
structure function, which tend to enhance the b induced process with respect to the g 
one. For example, for /x^ = m^, with rrit = 170(175) [180] GeV, one gets the 'leading 
logs' L = as(/i^) log(/i^/m^) ^ 0.75(0.76)[0.76]. Such differences between b and g ^ bb 



induced processes at the TeV scale are not unusual in literature, see, e.g., Refs. ^ 
(though, for the case of hadron-hadron collisions at the TeV scale). Note that we obtain 
the same pattern also for the case of on-shell top production, when no decay of the top 
quark is implemented. 

Before proceeding, we should in fact mention that we have studied the size of the 
differences between the total rates of the two processes as obtained, on the one hand, 
by using a finite width and implementing the decay currents and, on the other hand, by 
keeping the top on-shell. In general, they are at the level of few percent (the on-shell 
rates being larger). For example, for the b{g) induced process they vary between 2(1)% 
to 4(5)% when nit = 175 GeV. In fact, rates are rather insensitive to the value of the 
top mass. 

From Figure 2, we note in general that although the cross section for process (^) 
is small at existing collider energies {y/s^p ~ 300 GeV at DESY leads to a total cross 
section of less than 1 fb, which is negligible given the current integrated luminosity of 



about 20 pb ^ at each of the two experiments), it increases steeply near the TeV 



scale. At the LEP2©LHC scale it is easily observable at the 'conservative' luminosity of 



100 pb^ [^. There is however a sizeable dependence on the top quark mass, especially 
at low energies: the cross section being smaller for a phase space suppressed by a higher 
mass. (All our results hereafter assume the central value of 175 GeV.) 

Table I shows the cross section of the signal process (^ evaluated at the LEP2©LHC 
energy, for twenty-four different sets of structure functions. The dependence is found to 
be approximately 20%, with the maximum value of the total cross section differing from 
the minimum value by 823 fb. We believe such theoretical uncertainty to be already 
at the present time a reasonably small error so to motivate further and more detailed 
simulation studies (including hadronisation, detector effects, reducible background |]TB[) 
of single top phenomenology. To appreciate this we note that the result obtained by 
adopting the old set EHLQl |^ (i.e., the one used in Ref. |T6[) differs by that produced 
by MRS (A) in Table I by more than 50% ! 

As a further check, we present the Bjorken x and the Q dependence of the cross 
section of process @ for a selection of PDFs, in Figures 3a and 3b (respectively). In 
particular, we have included results for some older, MRS(Rl) and CTEQ(3M), and 
some newer, MRRS(l) and CTEQ(4M), sets, as representative of the two approaches 
MRS and CTEQ, each of the pair being fitted to a similar set of experimental data so 
to allow for a more consistent comparison. Note that the normalisations of the curves 
are to unity, in order to enlighten the differential behaviours of these quantities, in 
addition to their effects on the total rates (as was done in Table I). The clear message 
from Figures 3a and 3b is that the differences between the two pairs of sets are very 
small (as can be appreciated in detail in the central inserts), typically a few percent 
over all the available kinematic range in x and Q. Although we do not show the 
corresponding curves, we have verified that such considerations also apply to the other 
PDFs considered here. Thus, also at differential level the theoretical error on the rates 



of process (^) due to the PDFs is well under control already at present. 

Though it is beyond the scope of this study to trace back whether the differences 
in Table I (and Figures 3a-b) among the various sets are due to the gluon structure 
function or to the g ^ bb splitting (which onsets the b structure function), it is worth 
mentioning that it could well be that by the time a future ep collider will be running 
the uncertainties on the former will be so under control that one might attempt to 
distinguish between different dynamics proposed for the latter. In this respect, it would 
be interesting to assess whether the differences between MRRS(l) (dashed line) and 
CTEQ(4M) (dot-dashed line) in Table I and in Figures 3a and 3b are genuinely due 
to the dedicated treatment of the threshold region ~ /i;, performed in Ref. |3ll 



or 



not Clearly, this will require a tight control on all sources of experimental error, 
in particular of the actual value of the 6-tagging efficiency and of the hadronisation 
process of the quarks at the TeV scale. 

Figure 4 shows the differential distributions interesting for the final state phe- 
nomenology. Those in combined jet masses are sharply peaked at the top and 
masses well above the irreducible noise, indicating that for the hadronic case the jet 
masses can be used to clearly identify the top decays. This feature is convenient both 
for the elimination of top events from any other hadronic three-jet processes to which 
the former may act as a background, and for the elucidation of the top quark physics 
at ep colliders in, for example, probing the b quark distribution function. The spec- 
tra in transverse momenta Pt show that neither cuts in pt nor cuts in p™'^* affect 
the total cross section dramatically, whereas that of AR, the azimuthal-pseudorapidity 



separation defined by AR = y(A</))^ + (Arj)^ (where is the azimuthal angle and rj 
the pseudorapidity) indicates that the requirement of resolving the hadronic jets (or 
the requirement of an isolated lepton in the leptonic case) severely reduces the event 
rate. The majority of events are found within AR ~ 1.5, which is about 90 degrees in 
the azimuthal angle. This is because the visible jets and the lepton come from the en- 
ergetic top quark. Thus, at lower energies the azimuthal-pseudorapidity spread in the 
top quark decay products will be larger and hence the requirement of such jet/lepton 
isolation not so stringent. The distribution of the missing transverse momentum in the 
leptonic case, and more specifically the electronic case, is small at low missing pt and 
indicates that only a small proportion of the events will emulate neutral current events 
of the form ep eX. However, those which do will form a potentially dangerous 
background to high neutral current events, as the low AR mentioned above will 
concentrate the electrons to the high region. 

Table II and Figure 5 show the total cross section after the acceptance cuts. The fol- 
lowing LHC-like constraints were implemented (see |TB| for alternative selection strate- 



gies): for the leptonic channel, Pt ,Pt > 20 GeV, p^^"'^ > 10 GeV and AR£+^h > 0.7; 
for the hadronic channel, p^-' > 20 GeV, p™*''* > 10 GeV and ARjj'^h > 0.7. We 
have not introduced any cuts on pseudorapidity, as the particles were found to be all 
concentrated in the narrow \t]\ < 2.5 region even before any selection in pt was made. 

Table II summarises the event rates for all channels at three different CM energies. 
The numbers in square brackets are the cross sections of processes (P) and (|^). These 
are additional backgrounds when bottom quark charge tagging is unavailable. Since 
these effectively only differ from processes (^ and (^) in their non-resonant top quark 
production, they can be taken as a measure of the magnitude of the irreducible back- 



ground. As can be noticed, such background effects are small. We see that the hadronic 
cross section is higher at lower energy since, as discussed above, the acceptance cut in 
AR affects the rates less at smaller values of i/Sgp, thus compensating for the reduced 
total cross section shown in Figure 2. From Figure 5 we see that background effects 
do not spoil the sharp resonances in combined jet masses even after the acceptance 
cuts. We particularly stress that the QCD background is negligible: luckily enough, as 
it curiously peaks around the Mp(/± value in the di-jet mass distribution. It can also 
be noted that the cut in AR, the jet separation, of 0.7 is a conservative choice and, as 
can be deduced from the distribution in AR in Figure 4, so that the rates should be 
expected much higher for looser constraints, especially for the hadronic channel when 
three separate cuts in AR need to be made to resolve the three jets completely. Finally, 
given the optimistic vertex tagging performances foreseen for the LHC detectors 0, 
we would expect that only a small fraction of the event rates given in Table II will be 
lost in the actual analyses. In fact, we believe that the original capabilities of the LHC 
/i-vertex devices will be maintained while running the CERN machine in the proposed 
ep mode. 

4. Summary and conclusions 

The single top quark production from initial state bottom sea quark at future ep collid- 
ers was studied, mainly focusing our attention to the case of the proposed LEP2©LHC 
accelerator with the positron (electron) beam energy of 100 GeV and the proton one 
of 7 TeV. The total cross section was found to be about 4 pb at this energy. The 
uncertainty due to the structure functions was quantified to be rather small already at 
the present time, around 20%, and is expected to diminish significantly before new ep 
machines will enter into operation. Furthermore, based on such a consideration, some 
optimistic prospects about the possibility of exploiting single top phenomenology in 
order to study the g ^ bb dynamics inside the proton were given. Both the leptonic 
and hadronic decay channels of the top quark were studied, in presence of the corre- 
sponding irreducible backgrounds, which have been computed here for the first time. In 
the hadronic case, distributions in the reconstructed top quark and boson masses 
were found to be sharply peaked above the irreducible noise, so to allow for a prompt 
recognition of single top events. In both channels the cross section for the background 
was found to be small compared with the signal events. The residual dependence of the 
latter on the top quark mass was evaluated in several instances. Finally, the formula 
for the matrix element squared of single top production, including top width effects and 
all the dynamic correlations between the top decay products, was presented in order 
to aid future, more detailed, experimental simulations. The complete numerical pro- 
grams, evaluating irreducible background effects as well, have been especially optimised 
in view of high statistic Monte Carlo simulations and are available from the authors 
upon request. 
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Table Captions 



[I] Total cross sections (hadronic and leptonic) for process @ at LEP2©LHC energies 
for twenty-four different sets of structure functions. Errors are as given by VEGAS 
(the same statistics were used for the NCALL and ITMX parameters) p4| . 



[II] Total cross sections (hadronic and leptonic, including irreducible background ef- 
fects) for processes (^-(0) at the LEP2©LHC collider. The structure function set 



MRS (A) was used. Errors are as given by VEGAS p4|. The following acceptance 



cuts were implemented: (i) > 20 GeV, p^''' > 10 GeV and Ai?£+,fe > 0.7 

(leptonic channel); (ii) p'/,p^T > 20 GeV, p^''' > 10 GeV and ARjj,^b > 0.7 
(hadronic channel). In the squared brackets of the first two columns we report 
the rates of the charge conjugates (||)-(|^) of processes (^-(H), for which the res- 
onant top production do not occur. The rates of the charge conjugate of process 
(|^) are the same as those in third column. 



Figure Captions 



[1] Lowest order Feynman diagrams describing processes (|l|), (|), (H) and corre- 
sponding to sets (a), (b), (c) and (d), respectively. Only the cases l^vg, = e'^Ve 
and jj' = ud are shown are shown for processes whereas in reaction (|l|) 

the top is considered on-shell. The package MadGraph ||2l[] was used to produce 
the PostScript codes. In (c) 'A' represents the photon. The dashed lines in (b) 
and (c) represent the SM Higgs boson and the curly lines in (d) the gluon. The 
number of diagrams in (b) reduces to 10 for the cases £+1/^ = ji'^u^ and r"*"//,-, 
when diagrams 1, 2, 3, 5, 6, 7, 8, 12, 13, 20 and 21 do not contribute. 

[2] The total cross section (hadronic and leptonic channels) for processes (upper 
lines) and ([I|) (lower lines) for 300 GeV < v^ep — ^ TeV, with three different 
values for the top quark mass: = 170 GeV (continuous line), = 175 
GeV (dashed line) and = 180 GeV (dotted line). The structure function set 
MRS (A) was used. 

[3] Differential distributions in (a) x and (b) Q for events of the type (^) at the 
LEP2©LHC collider for three representative sets of structure functions: MRS(Rl) 
(solid), MRRS(l) (dashed), CTEQ(3M) (dotted) and CTEQ(4M) (dot-dashed). 
In the central inserts, the spectra are magnified around their maximum values. 
Normalisations are to unity. 

[4] Differential distributions (hadronic and leptonic channels) for process (^ at the 
LEP2©LHC CM energy and nit = 175 GeV in the following variables (clockwise). 
1. Mjcts, the invariant mass of the two- (solid) and three-jet (dashed) systems in 
hadronic decays. 2. px, the transverse momenta of the lepton/jets (solid) in 
leptonic/hadronic decays, of the bottom quark (dashed) in both channels, and 
of the missing particles in leptonic (dotted) and hadronic (dot-dashed) decays. 
3. AR, the azimuthal-pseudorapidity separation of the pairs lepton/jets-bottom 
quarks (solid) in leptonic/hadronic decays. The normalisation is to unity. The 
structure function set MRS (A) was used. In the case of the hadronic decays we 
have considered only one of the two light quark jets, their distributions in the 
above variables being very similar. 

[5] Differential distributions in Mj^t^ (hadronic channel only) for processes (^) (left) 
and (0) (right) at the LEP2©LHC CM energy and rrit = 175 GeV. M^,,, signifies 
the invariant mass of the two- (solid) and three-jet (dashed) systems in hadronic 
decays. The normalisations are to the total cross sections. The structure function 
set MRS (A) was used. Bins are 2 GeV wide. The following acceptance cuts were 
implemented: pj^-' ,p\^ > 20 GeV, p™'^'' > 10 GeV and ARjj/ h > 0.7 (hadronic 
channel) . 
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Table I 



o-tot (fb) 
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Diagrams by MadGraph 
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